sel1=c("M30_M29","M26_M25","M36_M35","M28_M27")
sel1=c("X2B_X1T","M8_M7","M6_M5","M2_M1","M42_M41","M44_M43","M48_M47","M50_M49","M52_M51")
sel1=c("M18_M17","M20_M19","M22_M21","M40_M39")
sel1=paste(sel1,".snp.ASM.sel.tsv",sep="")
countdiffa=c()
countdiffb=c()
countdiffc=c()
countdiffd=c()
for (i in 1:length(sel1)){
file=read.table(sel1[i],head=T,sep="\t")
count=paste(as.matrix(file)[,1],":",file[,2],sep="")
countdiffa=unique(c(countdiffa,count))
file=file[file$exits=="T",]
count=paste(as.matrix(file)[,1],":",file[,2],sep="")
countdiffb=unique(c(countdiffb,count))
if(dim(file)[1]!=0){
bias=(file$DU>=0.25 & file$p.biasU<0.05) | (file$DA>=0.25 & file$p.biasA<0.05)
diff=(file$p.diff<0.05)
countdiffc=unique(c(countdiffc,count[which(bias)]))
countdiffd=unique(c(countdiffd,count[which(bias & diff)]))
}
}